home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DPODI.z / DPODI
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPOOOODDDDIIII((((3333FFFF))))                                                            DDDDPPPPOOOODDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPODI   - DPODI computes the determinant and inverse of a certain double
  10.      precision symmetric positive definite matrix (see below) using the
  11.      factors computed by DPOCO, DPOFA or DQRDC.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE DPODI(A,LDA,N,DET,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAA DOUBLE PRECISION(LDA, N)
  21.         the output  A  from DPOCO or DPOFA
  22.         or the output  X  from DQRDC.
  23.  
  24.      LLLLDDDDAAAA INTEGER
  25.         the leading dimension of the array  A .
  26.  
  27.      NNNN INTEGER
  28.         the order of the matrix  A .
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         = 11   both determinant and inverse.
  32.         = 01   inverse only.
  33.         = 10   determinant only.  On Return
  34.  
  35.      AAAA If DPOCO or DPOFA was used to factor  A , then
  36.         DPODI produces the upper half of INVERSE(A) .
  37.         If DQRDC was used to decompose  X , then
  38.         DPODI produces the upper half of inverse(TRANS(X)*X)
  39.         where TRANS(X) is the transpose.
  40.         Elements of  A  below the diagonal are unchanged.
  41.         If the units digit of JOB is zero,  A  is unchanged.
  42.  
  43.      DDDDEEEETTTT DOUBLE PRECISION(2)
  44.         determinant of  A  or of  TRANS(X)*X  if requested.
  45.         Otherwise not referenced.
  46.         Determinant = DET(1) * 10.0**DET(2)
  47.         with  1.0 .LE. DET(1) .LT. 10.0
  48.         or  DET(1) .EQ. 0.0 .  Error Condition
  49.  
  50.      AAAA division by zero will occur if the input factor contains a zero on the
  51.      diagonal and the inverse is requested.  It will not occur if the
  52.      subroutines are called correctly and if DPOCO or DPOFA has set INFO .EQ.
  53.      0 .  LINPACK.  This version dated 08/14/78 .  Cleve Moler, University of
  54.      New Mexico, Argonne National Lab.  Subroutines and Functions BLAS
  55.      DAXPY,DSCAL Fortran MOD
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.